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I. INTRODUCTION 



The variational Monte Carlo (VMC) and Green's function Monte Carlo (GFMC) tech- 
niques are powerful tools for calculating properties of light nuclei. The GFMC method, 
in combination with the Argonne fig (AV18) two-nucleon (NN) and Illinois-2 (IL2) three- 



nucleon 



3iV) potentials, reproduces the ground- and excited-state energies for A < 10 nuclei 



[ll, 0, [3j_3]. It is beginning to be used to calculate reactions, such as nucleon- nucleus scat- 
tering 5|. Electroweak transitions in A = 6, 7 nuclei have been calculated using the more 
approximate VMC technique with AV18 and the older Urbana-IX (UIX) 3N potential. These 
include 6 Li elastic and transition form factors and radiative widths 0, electric quadrupole 
(E2) transition probabilities in 7 Li for pion inelastic scattering 0, Gamow-Teller (GT) ma- 
trix elements in 6 He and 7 Be weak decays |8|], and radiative capture reactions producing 6 Li, 
7 Li, and 7 Be @, M ■ 



The VMC results for electroweak transitions reported in earlier works were in reason- 
able agreement with experimental data. However, the AV18+IL2 Hamiltonian reproduces 
p-shell binding energies better than AV18+UIX, and GFMC wave functions are better ap- 
proximations to the true eigenstates. Hence GFMC calculations with AV18+IL2 for these 
electroweak transitions are worth investigating. In this work we study electromagnetic [E2 
and magnetic dipole (Ml)] transition strengths and nuclear beta-decay [Fermi (F) and GT] 
rates, with the GFMC technique using the AV18+IL2 potential for nuclei with A = 6, 7. 
This is the first off-diagonal matrix element calculation using the nuclear GFMC method. 
The GFMC technique of evaluating off-diagonal matrix elements should be applicable in 
many other nuclear calculations, such as isospin-mixing, and low-energy nuclear reactions. 

In this work we consider only the one-body parts of the transition operators. Schematic 
expressions for these are 



E2 = ^[rM](l + r b ), (1) 

k 

Ml = fi N J2[(Lk + 9 P S k )(l + r kz )/2 + g n S k (l - r kz )/2], (2) 

k 

F = E^±> (3) 



k 



GT = 5> fc r fc± , (4) 

k 

where k labels individual nucleons, r is the spatial coordinate, Y is a spherical harmonic, 
7fez(±) is the third (raising/lowering) component of the isospin operator, a is the Pauli spin 
matrix, fi^ is the nuclear magneton, L(S) is the orbital (spin) angular momentum operator, 
and g p ( n ) is the gyromagnetic ratio for a proton (neutron). In all cases the appropriate 
projection (z, +, and — ) of the operators is understood. Note that e, fi^, and g p ( n ) are all 
physical values, i.e., no effective charge is used. 

The current work does not include the contributions of two-body electroweak current op- 
erators. The effect of two-body meson-exchange operators has been studied previously, and 
found to be small in electric quadrupole and weak transitions, but significant for magnetic 
transitions (1, 
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II. VMC TRIAL FUNCTIONS 



In this work we calculate the off-diagonal transition matrix element 
{^/f(J 7T/ ; T")|0|\E f j(J 7r ; T)), where O is one of the one-body electroweak transition op- 
erators given above and \&( J n ; T) is the nuclear wave function with a specific spin-parity J n 
and isospin T. The variational wave function is an approximate solution of the many-body 
Schrodinger equation 



Hty{J w ;T) = E^{J n ;T) . 
The Hamiltonian used here has the form 



(5) 



'j 



i<j 



i<j<k 



(6) 



with Ki being the non-relativistic kinetic energy, Argonne v\& [H| is the NN potential v 



and Illinois-2 12 



is the 3N interaction V^. 
The VMC trial function ^T{J n ', T) for a given nucleus is constructed from products of 
two- and three-body correlation operators acting on an antisymmetric single-particle state 
with the appropriate quantum numbers. The correlation operators are designed to reflect the 
influence of the interactions at short distances, while appropriate boundary conditions are 
imposed at long range 0, [Io|, 13, 14]. The ^(J^T) has embedded variational parameters 
that are adjusted to minimize the expectation value 



(\f! T \\f/ T ) 



• 



(7) 



which is evaluated by Metropolis Monte Carlo integration 15 
A good variational trial function has the form 



|*t) 



i<j<k 



TNI 

ijk 



511(1 



(8) 



The Jastrow wave function, tyj, is fully antisymmetric and has the (J n ;T) quantum numbers 
of the state of interest. For s-shell nuclei, \Pj has the simple form 



11 fijk 




n fMj) 


i<j<k<A 




i<j<A 



\<$> A {JMjTT z )) 



(9) 



Here f c {T%j) and f^ k are central two- and three-body correlation functions and &a is a Slater 
determinant in spin-isospin space, e.g., for the a-particle, |$ 4 (0000)) = A\p t P I n t n I}- 
The Uij and U™ 1 are noncommuting two- and three-nucleon correlation operators, and 
S indicates a symmetric product over all possible ordering of pairs and triples. The 
includes spin, isospin, and tensor terms: 



ij 



p=2,6 



(10) 



where the O 



P=l,6 
ij 



[l,(jj • (Tj,Sij) ® [l,Tj ■ Tj] are the main static operators that appear 
in the NN potential. The f c (f) and u p (r) functions are generated by the solution of a 
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set of coupled differential equations containing the bare NN potential with asymptotically- 



confined boundary conditions [13j . The Uf^ 1 has the spin-isospin structure of the dominant 



ijk 

parts of the 3N interaction as suggested by perturbation theory. 

For the p-shell nuclei, *S?j includes a one-body part that consists of 4 nucleons in an a-like 
core and (A — 4) nucleons in p-shell orbitals. We use LS coupling to obtain the desired JMj 
value, as suggested by standard shell-model studies [l6| . We also need to sum over different 
spatial symmetries [n] of the angular momentum coupling of the p-shell nucleons [17J ]. The 
one-body parts are multiplied by products of central pair and triplet correlation functions, 
which depend upon the shells (s or p) occupied by the particles and on the LS[n] coupling: 

i*j> = m n fw n /•.(»•«) n 

i<j<k i<j<4 k<4<l<A 

E (Plsw II / p L / W (n m )|$A( J L5HJM J TT 3 ) 12 34:5...A)) } • (11) 

LS[n] 4<Km<A 

The operator A indicates an antisymmetric sum over all possible partitions into 4 s-shell 
and (A — 4) p-shell particles. 

The LS'jn] components of the single-particle wave function are given by: 



\® A (LS[n}JMjTT 3 ) 12M:5 ... A ) = |$ 4 (0000) 1234 II 0f W (^ 

4<1<A 



1 1 

II Yim^cd^LMM x [ II Xi(~m s )]s Ms \ x [ H vi{-h)) TTa ) ■ (12) 

4<1<A 4<1<A Z ) JM 4<Z<A Z 

The <f)p S (R a i) are p-wave solutions of a particle in an effective a-N potential that has Woods- 
Saxon and Coulomb parts. They are functions of the distance between the center of mass of 
the a core and nucleon I, and may vary with LS[n}. The f ss , f sp , and fp P s ^ all have similar 
short-range behavior, like the f c of the a-particle, but different long-range tails. 

Two different types of have been constructed in recent VMC calculations of light 



p-shell nuclei: an original shell-model kind of trial function [14j which we will call Type I, 
and a cluster-cluster kind of trial function 0, 10] which we will call Type II. In Type I trial 



functions, the (fip S (r) has an exponential decay at long range, with the depth, range, and 
surface thickness of the Woods-Saxon potential serving as variational parameters. The f sp 
goes to a constant ~ 1, while fp p s ^ has a much smaller tail to allow clusterization of the 
p-shell nucleons. Details of these A = 6, 7 trial functions are given in Ref. [14J . 

In Type II trial functions, 4>p S {r) is again the solution of a p-wave differential equation 
with a potential containing Woods-Saxon and Coulomb terms, but with an added Lagrange 
multiplier that turns on at long range. This Lagrange multiplier imposes the boundary 
condition 

[<p L p S[n] (r - oc)p- 4 ) oc ^ fcm (2 7 r)/r, (13) 

where the Whittaker function, Wkmi^jr), gives the asymptotic form of the bound-state wave 
function in a Coulomb potential. The 7 is related to the cluster separation energy which is 
taken from experiment (the GFMC computed separation energies for AV18+IL2 are close 
to the experimental values). The accompanying f sp goes to unity (more rapidly than in the 
Type I trial function) and the f pp ^ are taken from the exact deuteron wave function in the 
case of 6 Li, or the VMC triton ( 3 He) trial function in the case of 7 Li ( 7 Be). Consequently, 
the Type II trial function factorizes at large cluster separations as 

-> ipa^TW km (2^r aT )/r aT . (14) 
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where ip a and ijj T are the wave functions of the clusters and r aT is the separation between 
them. More details on these wave functions are given in Refs. 0, [Toj]. In the case of 6 He, 
which does not have an asymptotic two-cluster threshold, we generate a /^f^ correlation 
assuming a weakly-bound 1 S nn pair. 

For either type of trial function, a diagonalization is carried out in the one-body basis to 
find the optimal values of the ^LS[n] mixing parameters for a given (J n ;T) state. The trial 
function, Eq.©, has the advantage of being efficient to evaluate while including the bulk of 
the correlation effects. 



III. GFMC WAVE FUNCTIONS 



The GFMC method 18, 19| projects out the exact lowest-energy state, for a given 
set of quantum numbers, using 



*o = lim exp[-(H' - E )T]^q 



(15) 



where H' is a possibly simplified version of the desired Hamiltonian H and \I/t is an initial 
trial wave function. If the maximum r actually used is large enough, the eigenvalue Eq is 
calculated exactly while other expectation values are generally calculated neglecting terms 



of order 
of order 



^t\ 2 and higher [14|. In contrast, the error in the variational energy Ey is 



^t\ 2 , and other expectation values calculated with ty? have errors of order 



\jr _ \]> r | j n th e following we present a brief overview of nuclear GFMC methods; much 



more detail may be found in Refs. [14|, [20 . 

We start with the of Eq.([H]) and define the propagated wave function \I/(t) 



-{H'~E )At 



(16) 



where we have introduced a small time step, r = nAr; obviously \l/(r = 0) = and 
\I/(t — > oo) = ^o- Quantities of interest are evaluated in terms of a "mixed" expectation 
value between and \&(r): 



(0(t)) 



M 



(*(r)|0|* T > 



(17) 



(*(r)|* r > ' 

The desired expectation values would, of course, have ^/(r) on both sides; by writing ^(r) = 
+ S^(t) and neglecting terms of order [<5\I/(t)] 2 , we obtain the approximate expression 



(O(t)) 



(*(r)|0|*(r)> 



(0(t)) m + [{0(t)) m - (0) v ] 



(18) 



(*(r)|*(r)> 

where ( 0)v is the variational expectation value. More accurate evaluations of (0(t)) are 
possible [2l| , essentially by measuring the observable at the mid-point of the GFMC propa- 
gation. However, such estimates require a propagation twice as long as the mixed estimate 
and require separate propagations for every expectation value to be evaluated. The nuclear 
calculations published to date use the approximation of Eq. ffTSl . 

For off-diagonal matrix elements relevant to this work the mixed estimate is generalized 
to the following expression 

(^ / (r)IO|^ i (r)) 

' ' ' ! 1 " - (0(r)) Mi + (0{r)) Mf - (0) v , (19) 
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where 



(Oh 



[^\0\¥ T ) 



(0(t)) 



M, 





0|**(r)) 


(*VI^> 




P(r)> \ 







(0(r)> 



Mr 



(^f(r)\0\^ T ) 

mr)\H) \ 



(20) 



(21) 



(22) 



and the index i (f) refers to the wave function of the initial (final) nuclear state. In our 
calculation the operator always acts on the trial wave function. The first term of Eq. (|2ip 
is thus replaced by its conjug ate <**(T)|Ot|^/) 

in our real computation. Note here the 
computation for each extrapolated transition requires two independent GFMC propagations. 

The quantities (^(t)|0|^)/(^(t)|^) and (t)\0\W t ) / (r)\ty f T ) can be directly 
evaluated in the GFMC propagations of \l/ l (r) and ^(r) respectively. The VMC expectation 
value can be cast as 



(O) 



v 



) \ 



T 



o\H) 



(nl^r) 



(23) 



(24) 



in which the first term can be computed in VMC walks guided by ty l T and respectively. 
The ratio (*tI*t>/(*tI*t) or its 

inverse can be also computed in the ^/ l T and walks. 
In this present GFMC calculation the propagation of the mixed off-diagonal matrix element 
was generally carried out for a value of r up to 3 MeV -1 . In most cases the transition matrix 
elements are quite stable with this large value of r. 

As noted in Eq. ffl5|) . the GFMC propagation may be computed using a simplified Hamil- 
tonian. In the current work we use H' = AV8' + IL2' where AV8' is a reprojection of AV18 
defined in Ref. [13] and the strength of the central repulsive part of IL2 is modified in IL2' so 
that (if') w (H). Energies are then perturbatively corrected by adding (H — H'). However 
all other expectation values are really for eigenfunctions of H'\ we have no way of correcting 
them to expectation values in the eigenfunctions of the desired H. Because H' is a good ap- 
proximation to H, this in general should not be a problem. However, Fermi matrix elements 
are different from their trivial values (2J+1), only because of charge-independence-breaking 
components in the wave functions: 



(2J+l)(l-e) 



(25) 



The AV8' does not contain the strong charge symmetry breaking (CSB) component of AV18 
and has only a Z-dependent isoscalar projection of the Coulomb potential as proposed by 
Kamuntavicius, et al. (22|; all other electromagnetic terms in AV18 are not included. There- 
fore our calculations may seriously underestimate the correct values of e for AV18+IL2. This 
is suggested by a comparison with correlated hyperspherical harmonics (CHH) calculations, 
using the AV18+UIX Hamiltonian, of the F and GT matrix elements for 3 H decay j23|. The 
CHH value for the GT matrix element is 2.258 while our result computed using AVS'+UIX' 
is 2.260±0.001. However their value for e is 0.0013, while ours is 0. ±0.0005. 
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TABLE I: Electromagnetic transition reduced matrix elements for A = 6 nuclei. 



,J l ! X l 


>JJ;T f 


mode 


VMC 


Mi 


M f 


GFMC 


6 Li(3+;0) - 
6 Li(2+;0) - 
6 Li(0+;l) - 
6 Li(2+;l)- 


^ 6 Li(l+;0) 
^ 6 Li(l+;0) 
^ 6 Li(l+;0) 
^ 6 Li(l+;0) 


E2(e fm 2 ) 
E2(e fm 2 ) 
Ml(n N ) 
Ml(n N ) 


8.20(1) 
6.20(8) 
3.682(4) 
0.09(3) 


8.46(3) 
6.49(16) 
3.658(1) 
-0.010(3) 


8.77(3) 
6.29(3) 
3.643(1) 
-0.010(7) 


9.03(5) 
6.58(18) 
3.619(4) 
-0.11(3) 


6 He(2+;l) - 


^ 6 He(0+;l) 


E2{e fm 2 ) 


1.44(1) 


1.64(2) 


1.63(1) 


1.81(3) 



IV. RESULTS FOR A = 6, 7 NUCLEI 

Evaluation of the transition matrix elements is fairly straightforward. The VMC or 
GFMC wave function samples for a given (J 77 ; T) state are constructed with a specific Mj 
projection. We find it convenient to use the same Mj for both initial and final states, even 
if J is different. This makes the Ml operator exactly equivalent to the magnetic moment 
operator, and also keeps the GT operator particularly simple. We note that the size of the 
Monte Carlo statistical errors can vary significantly depending on the particular Mj substate 
that is chosen. Below we present more than a dozen electroweak transitions between different 
states of A = 6, 7 nuclei. The first two subsections discuss the electromagnetic transitions 
and the last subsection discusses the weak transitions. 



A. Electromagnetic Transitions of A = 6 Nuclei 

Table [JJ shows the VMC, the two mixed estimates, and the extrapolated GFMC re- 
duced matrix elements given in Eq. (jl9p . for various electromagnetic transitions between 
different states of A = 6 nuclei. The corresponding transition widths (computed with the 
experimental excitation energies) are shown in Table [Til where they are compared with 
Cohen-Kurath (CK) shell-model values [l6j], no-core shell model values (NCSM) [241], and 



experiment [25|, |26j. The NCSM values were computed for the AV8'+TM' Hamiltonian 
which we expect to have similar transitions to the AV8'+IL2' Hamiltonian used here. 

The GFMC propagation for the 6 Li E2 matrix elements is shown in Fig. [TJ For each of 
the E2 transitions we plot the two mixed reduced matrix elements: red squares for those 
with GFMC ground state configurations and green circles for those with GFMC excited state 
configurations. The solid purple line starting from the origin shows the pure VMC estimate, 
while the black stars represent the extrapolated matrix elements in each E2 transition. The 
other solid lines are the average, over the range of r shown, for each reduced matrix element, 
with standard deviations shown as dashed lines. 

In Table[TJ and Fig{T]we note that for the E2 transition between the ground (1 + ; 0) and 
the first excited (3 + ;0) states of 6 Li, the average values for the two mixed estimates are 
larger than the pure VMC estimate. As a result the extrapolated GFMC matrix element is 
10% larger than the VMC value and hence the GFMC transition width is about 20% larger 
than the VMC width, as shown in Table [UJ We also note that the GFMC width for this 
transition is a little bigger than the experimental value, but it is within the experimental 
range. It is worth mentioning that the CK shell-model prediction for this width is about 
half the experimental value, despite the use of effective charges for proton and neutron of 



7 



TABLE II: Electromagnetic transition widths in eV of A = 6 nuclei. 



,J l ! X l 


>JJ;T f 


EM mode 


CK 


NCSM a 


VMC 


GFMC 


Expt. 


6 Li(3+;0) - 


^ 6 Li(l+;0) 


El (10" 4 ) 


2.17 


1.22 


3.86(1) 


4.68(5) 


4.40(34) 


6 Li(2+;0) - 


^ 6 Li(l+;0) 


El (10~ 2 ) 


0.34 


0.38 


0.92(2) 


1.04(6) 


0.54(28) 


6 Li(0+;l) - 


^ 6 Li(l+;0) 


Ml (10°) 


7.85 


8.05 


7.10(2) 


6.86(2) 


8.19(17) 


6 Li(2+;l)- 


^ 6 Li(l+;0) 


Ml (10°) 


0.23 


0.20 


0.003(2) 


0.004(2) 


0.27(5) 


6 He(2+;l) - 


^ 6 He(0+;l) 


El (1(T 5 ) 






0.63(1) 


1.03(3) 


1.63(10) 



a for AV8'+TM' from Ref. [24 1 
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FIG. 1: (Color online) E2 transitions for 6 Li(3+;0) and 6 Li(2 + ;0) to 6 Li(l+;0) ground state. 



1.4e and 0.4e, respectively. The NCSM value is only one quarter our result, despite our 
expectation that the two Hamiltonians (AV8'+TM' in NCSM and AV8'+IL2' in our work) 
should not result in very different transition moments. In a more recent publication, NCSM 
values using a IQTiu space for AV8' with no 3N potential were presented [27j|. The B(E2) 
increased by 20% from the 6hu values presented in Table fl] 

The E2 transition from (2 + ; 0) excited state to the ground state turns out to be difficult 
to calculate, as can be seen in Fig. [IJ The plot shows that the GFMC mixed estimates using 
either the ground or the first excited (3 + ; 0) states of 6 Li nucleus are moderately stable. 
On the other hand the GFMC mixed estimate which uses the broad (2 + ; 0) state of 6 Li is 
growing rapidly with r, which makes a simple average meaningless. This is undoubtedly 
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0.5 

FIG. 2: (Color online) GFMC E(r) for different states of 6 Li. 

related to the difficulty in obtaining GFMC energies for broad, particle-unstable, states. 
Figure [2] shows the computed energies of the (l + ;0), (3 + ;0), and (2 + ;0) states in 6 Li as a 
function of r. The (l + ;0) and (3 + ; 0) energies drop rapidly with r from the initial VMC 
value and then become constant, aside from statistical fluctuations. The stable energy is 
reached around r = 0.2 MeV -1 as marked in the figure by the open stars. However, after a 
similar initial rapid decrease, the energy of the experimentally broad (2 + ; 0) state continues 
to decrease. At the same time the rms radius is steadily increasing; the GFMC algorithm 
is propagating this state to separated a and deuteron clusters. Based on the convergence 
of the (l + ;0) and (3 + ;0) energies, we assume that values at r = 0.2 MeV -1 represent the 
best GFMC estimates for this state. Using this value of r for the transition results in the 
value 6.49(16) for the (2 + | \E2\ |1 + ) Mi matrix element, which is shown as an open circle 
in Figure [TJ The quoted error is based on a range of ±0.1 for the r at which the value is 
evaluated. This GFMC result for the width is only 10% larger than the VMC value, but 
is three times as big as the CK value, and twice as big as the experimental value, which 
however has a sizeable error bar. 

A similar analysis was used for the 6 He E2 matrix element given in Table [B The experi- 
mental value here is taken from a recent measurement of the B(E2 from 6 He breakup on 
209gj near the Coulomb barrier 26| . 

The GFMC propagation for the Ml transition 6 Li(0 + ; 1) — > 6 Li(l + ; 0) is shown in Fig. 
The GFMC matrix element reduces the VMC estimate slightly, giving a width that is smaller 
than the current experimental value. However, this is not unexpected as we have used 
only one-body transition operators in our present calculation. Two-body meson-exchange 
currents (MEC) are known to increase isovector magnetic moments by 15-20% for A = 3 
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FIG. 3: (Color online) Ml transition for 6 Li(0+; 1) to 6 Li(l+;0). 

nuclei, while having profound effects on the magnetic form factors [28]. A previous VMC 
calculation of the width of this transition found a 20% increase from 7.49(2) eV to 9.06(7) 
eV when MEC appropriate for the AV18+UIX Hamiltonian were added [6j. A similar 
increase applied to the present GFMC calculation would predict a width of 8.29(3) eV in 
excellent agreement with the experimental width of 8.19(17) eV. We plan to evaluate MEC 
corrections with the GFMC wave functions in future work. The NCSM result is in good 
agreement with experiment without any MEC contributions; if the MEC contributions are 
the size we expect, however, this good agreement will be lost when they are added. 

Finally, we have a very difficult time evaluating the Ml transition between 6 Li(2 + ; 1) 
and 6 Li(l + ; 0) states. The former is again a wide state which has the same GFMC propaga- 
tion difficulty as the E2 transition between 6 Li(2 + ; 0) and 6 Li(l + ; 0) states discussed above. 
However, the biggest problem is a large cancellation between different LS[n] components of 
the two wave functions, because the dominant pieces, 1 D 2 [42] and 3 Si[42], are not connected 
by the Ml operator. The VMC diagonalization and GFMC propagation are driven to mini- 
mize the energy, and may not determine small components of the wave functions sufficiently 
well to obtain such sensitive cancellations. Finally, the contribution of MEC terms may be 
much more important here because of the cancellations in the impulse approximation; hints 
of this were observed in the earlier VMC study 
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TABLE III: Electromagnetic transition reduced matrix elements of A = 7 nuclei. 



J 


i 






mode 


VMC 


M t 


M f 


GFMC 


7 Li(|" 


> 


^Li(f-) 


I 


E2(e fm 2 ) 


5.11(5) 


5.44(2) 


5.37(2) 


5.69(6) 


?Li(I- 


")- 


^Li(l-) 


II 


E2(e fm 2 ) 


5.38(6) 


5.53(3) 


5.56(2) 


5.71(7) 




") " 


^ 7 Li(§ ) 


I 


Ml(/ijv) 


2.742(1) 


2.749(3) 


2.693(2) 


2.695(4) 


7 Li(I" 


")" 


^Li(r) 


II 


Ml(/ijv) 


2.738(1) 


2.673(3) 


2.706(2) 


2.641(3) 




')" 


- 7 Li(l") 


I 


E2(e fm 2 ) 


7.67(4) 


8.28(3) 


8.30(3) 


8.91(6) 


7 Be(i 


")- 


- 7 Be(l") 


I 


£2(e fm 2 ) 


8.51(3) 


9.09(4) 


9.54(3) 


10.12(6) 


7 Be(i" 


')" 


- 7 Be(f-) 


II 


£2(e fm 2 ) 


8.86(13) 


9.74(7) 


9.60(6) 


10.48(16) 


7 Be(i" 


')" 


- 7 Be(f-) 


I 


Ml(p N ) 


2.423(2) 


2.412(2) 


2.403(3) 


2.394(4) 


7 Be(i" 


> 


- 7 Be(l") 


II 


Ml(ji N ) 


2.405(3) 


2.390(2) 


2.386(5) 


2.372(6) 



B. Electromagnetic Transitions of A = 7 Nuclei 

In Table ILTT1 we present the matrix elements of a number of electromagnetic transitions in 
A = 7 nuclei. As in Table [U this table shows the VMC estimates, the two mixed estimates 
and the GFMC extrapolated matrix elements for each transition. We suppress the isospin 
quantum numbers for different states of 7 Li and 7 Be because all states we consider have 
T — |. For those transitions in A = 7 nuclei between particle-stable states, we made two 
independent calculations using both Type I and Type II trial wave functions as discussed 
in Sec. [TTJ However, the Type II trial function is not defined for particle-unstable states like 
7 Li(| ). It is expected that, even though the VMC estimates may be somewhat dependent 
on the trial wave functions, the GFMC calculation should remove most of the dependence. 
(It is exact at the order of j^o — ^t| 2 )- The extrapolated regular (I) and asymptotic (II) 
expectation values are within 2% of each other for every transition we considered. 

Table [IV] shows the corresponding widths of the electromagnetic transitions in A = 7 
nuclei, compared to the CK shell-model values 16] and experiment 25, 2^]. The GFMC 



width is about 20% bigger than the CK and VMC values, and in good agreement with the 

:f)- 7 Li(§- 



experimental width for the 7 Li(^ ) ^ 7 Li(| ) E2 transition. The corresponding transition 



in 7 Be has not been measured. The widths for both the Ml transitions are relatively smaller 
than the corresponding experimental widths as expected for just one-body magnetic-moment 
operator expectation values. If there is a 20% additional contribution as expected from MEC 
terms, these will approach within 10% of the experimental values. 

The GFMC propagation for the E2 transition from 7 Li(| ) to 7 Li(| ) is illustrated in 
Figure HI We note that even though the individual points for the two mixed estimates are 
a little scattered, the average values for these overlap. The extrapolated GFMC result is 
larger than the VMC estimate by 15%, making the transition width 30% larger. The VMC 
value is one experimental standard deviation below the experimental value 29[ while the 
GFMC value is the same amount above the experimental value. 
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TABLE IV: Electromagnetic transition widths in eV of A = 7 nuclei. 



J7 -» J? 



EM mode 



CK 



VMC 



GFMC 



Expt. 



7 Li(I-) - 


^Li(f-) 


I 


E2 (1(T 7 ) 


2.79 


2.61(3) 


3.24(7) 


3.30(20) 


? Li(D " 


^Li(l-) 


II 


£2 (icr 7 ) 




2.90(3) 


3.26(8) 


3.30(20) 


->'Li(§") 
^Li(f-) 


I 


Ml (W- 3 ) 


5.69 


4.74(3) 


4.58(3) 


6.30(31) 


? Li(D - 


II 


Ml (10~ 3 ) 




4.73(1) 


4.41(1) 


6.30(31) 


? Li(I") " 




I 


E2 (1(T 2 ) 


0.98 


1.29(1) 


1.74(2) 


1.50(20) 



7 Be(H 


- 7 Be(l-) 


I 


E2 (1Q- 7 ) 




4.24(3) 


6.00(7) 




7 Be(H 


- 7 Be(§') 


II 


E2 (10- 7 ) 




4.60(13) 


6.44(19) 




7 Be(H 


- 7 Be(f-) 


I 


Ml (10~ 3 ) 




2.69(1) 


2.62(1) 


3.43(45) 


7 Be(H 


- 7 Be(l") 


II 


Ml (10" 3 ) 




2.65(1) 


2.57(1) 


3.43(45) 





11.0 
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10.0 
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eg 




9.0 








8.5 






a 


8.0 


^> 




V 




7.5 
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FIG. 4: (Color online) E2 transition for 7 Li(| ) to 7 Li(|" 



C. Weak Transitions in A = 6, 7 nuclei 

The weak Fermi and Gamow- Teller matrix elements in A = 6, 7 nuclei are shown in Ta- 
ble [VJ We note that in every F and GT transition the extrapolated GFMC matrix elements 
are smaller than the VMC estimates. However, the reduction is not large, being about 2% 
for GT terms. This suggests the starting trial functions are already good approximations 
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TABLE V: Weak transition reduced matrix elements of A = 6, 7 nuclei. 



T7T 

J i 




\T/ , , 




V 1V-LV^ 


IVl t 


Mi 

IVlf 


V.J 1 IV J. 


nci u , ± i 




C\\ T 








9 1 7^1 \ 


9 1 1 ^ 


6 He(0 + ; 1) 


^ 6 Li(l+ 


0) II 


GT 


2.253(3) 


2.184(1) 


2.276(1) 


2.207(3) 


7 Be(|") 


- 7 Li(i" 


) I 


F 


2.0000(0) 


1.9998(1) 


1.9998(1) 


1.9997(3) 


7 Be(|") 


- 7 Li(F 


) II 


F 


1.9995(1) 


1.9987(3) 


1.9983(3) 


1.9976(5) 


7 Be(|") 


^ 7 Li(|- 


) I 


GT 


2.325(1) 


2.298(1) 


2.301(1) 


2.274(2) 


7 Be(|") 


^Li(|- 


) II 


GT 


2.339(4) 


2.311(1) 


2.319(1) 


2.291(4) 


7 Be(|") 


^Li(I- 


) I 


GT 


2.146(2) 


2.119(3) 


2.129(3) 


2.099(4) 


7 Be(|") 


- 7 Li(F 


) II 


GT 


2.139(1) 


2.121(1) 


2.098(2) 


2.080(1) 


7 Be(±") 


- 7 Li(F 


) II 


GT 


2.138(1) 


2.125(3) 


2.104(1) 


2.092(3) 



for these weak transitions. The differences between Type I and Type II trial functions are 
not great, and the GFMC propagation does not reduce these differences. 

The Fermi matrix element for 7 Be(| ) to 7 Li(| ) is exactly 2 for charge-symmetric wave 
functions such as our Type I. The Type II wave functions are not charge-symmetric because 
the at and ar separation energies in Eq. (fl~4l are different and also the triton and 3 He clusters 
are slightly different. It appears that the GFMC propagation will introduce only a small 
asymmetry when starting from a charge-symmetric (Type I) trial function, but if given a 
small starting asymmetry (Type II), it can enhance it considerably. However, as noted at 
the end of Sec. Ill, the present calculations may still seriously underestimate this asymmetry 
and its effect on We also show in the last line of Table IVl the GT matrix element for 

the transition from 7 Be(| ) to 7 Li(| ) which is not an observable weak decay, but could be 
measured in a (p, n) reaction. The charge-symmetric Type I trial function would give the 
same result as the 7 Be(| ) to 7 Li(| ) transition, but the Type II trial function gives a very 
slightly different result. 

The log(/t) values obtained from VMC, shown in Table IVIl are already in reasonable 
agreement with the corresponding experimental values and the GFMC values are even better. 
The previous VMC study included MEC contributions which boosted the GT transition 
matrix elements ~ 1.5% for A = 6, and ~ 3% for A = 7. This resulted in too small a 
half-life for 6 He but about right for 7 Be. When MEC contributions are eventually added to 
the GFMC calculation, the half-life for 6 He should be quite good, but the rate for 7 Be will 
probably be a little too fast. The last two lines of Table [VI] give the branching ratio £ of 
the weak decay to the two final states in 7 Li for Type I and II trial functions. These are 
also a little low compared to experiment, but MEC contributions should also improve the 
agreement. 

In Fig. we present the various reduced matrix elements as a function of r for two 
Gamow- Teller transitions, 7 Be(| ) to 7 Li(| ) and 7 Be(| ) to 7 Li(| ). The former is shown 
for the Type I trial function, and the latter for the Type II. The GFMC mixed estimate 
points for both transitions are quite stable. 
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TABLE VI: Log(/t) values for weak transitions of A = 6, 7 nuclei 



T7T 
J i 


, 77T 




Weak Current 




in Uoivi 


V 1V1L/ 




i ,' . . , , * 

HjXpi. 


ne^u , i ) 


— > L/l^l , U) 


T 
1 


prp 


9 Si A 


9 87 
Z.ol 


9 Q01 ^ 




9 Q1 0^9^ 


6 He(0+;l) 


-» 6 Li(l+;0) 


II 


GT 






2.879(2) 


2.897(2) 




7 Be(|-) 


^Li(f-) 


I 


F & GT 


3.38 


3.30 


3.288(1) 


3.302(1) 


3.32 


7 Be(i~) 


- ? Li(l") 


II 


F & GT 






3.285(1) 


3.297(1) 




7 Be(D 




I 


GT 


3.46 


3.53 


3.523(1) 


3.542(1) 


3.55 


7 Be(|-) 


- 7 Li(D 


II 


GT 






3.526(1) 


3.550(1) 




? Li(i" 


)M1 ) 


I 


F & GT 


14.2% 


10.38% 


10.38(3)% 


10.25(3)% 


10.44(4)% 


Z Li(I" 


)/Li(|~) 


II 


F & GT 






10.25(3)% 


10.00(3)% 





Q for AV8'+TM' from Ref. [2|] 



2.50 



2.40 - 



~ i 1 1 1 1 1 1 1 — i — | — i 1 1 1 1 1 1 1 1 — | — i — i 1 1 1 1 1 1 r 

■ Mixed f 7 7. - 2.50 
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3 ,... j. ... ; ... A ... S .... T .... $ .... ft ...; 



I 



5 

3/2" -> 1/2 




2.30 



- 2.20 



T I 



2.10 -S 



2.10 



2.00 
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I 



2.00 
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x (MeV ) 



FIG. 5: (Color online) Gamow-Teller transitions for 7 Be(§ ) to 7 Li(§ ) and 7 Be(§ ) to 7 Li(± ). 



V. CONCLUSIONS 

These first GFMC calculations of transition matrix elements in light A = 6, 7 nuclei 
are generally in good agreement with the current experimental data. A number of these 
transitions have been calculated previously using the more approximate VMC technique 
with an older 3 A" potential 0, 0, H| ■ Here we explored a significant number of electroweak 
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transitions using the GFMC method, and calculated the corresponding widths or log(/£) for 
each transition. We compared our results to Cohen-Kurath shell-model and no-core shell- 
model results where available, and also with the current experimental results. In most of the 
transitions we considered we found that the GFMC transition widths or log(/i) values have 
been improved from the VMC and are in good agreement with experimental numbers. This 
is in general true for most of the E2 and all the GT and F type transitions. However, for Ml 
type transitions the GFMC widths we obtained are smaller than the current experimental 
values. Meson-exchange current corrections are expected to be large for Ml transitions and 
must be calculated for a meaningful comparison with data. We note here that the effect 
of MEC on E2 and GT transitions are expected to be smaller, about < 3%. In addition 
to the good results we obtained in most cases we faced some difficulties, especially treating 
broad nuclear states using the GFMC method. In these cases scattering boundary conditions 
should be used; GFMC has recently been successfully applied to the na scattering states jl] . 
We also had difficulty when the main components of the wave functions did not contribute 
to the transition, with the result depending on cancellations between small components. 

Some of the transitions that we explored were also treated by using two different types of 
trial wave functions. The extrapolated GFMC values obtained by using either of the wave 
functions should be the same; in practice they are within 2% of each other. This is indeed 
found in our calculations; the widths we obtain by using one or the other type of trial wave 
function are very close. 

In future, we expect to extend this work to larger nuclei in the ^4=8-10 range and for 
additional operators such as El and M3. One difficulty we anticipate is that some tran- 
sitions of interest, such as the weak decays of 8 He, 8 Li, and 8 B, run predominantly from 
large components of the initial state wave function to small components in the final states. 
These small components may not be well-determined by the GFMC calculation so addi- 
tional constraints may be necessary. We also need to evaluate two-body contributions to 
the electroweak current operators that are consistent with our chosen Hamiltonian. 
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